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We propose a new algorithm MARINLINGA for reverse line graph computation, i.e., construct- 
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ing the original graph from a given line graph. Based on the completely new and simpler principle 
of link relabeling and endnode recognition, MARINLINGA does not rely on Whitney's theorem 
while all previous algorithms do. MARINLINGA has a worst case complexity of 0{N^), where 
iV" denotes the number of nodes of the line graph. We demonstrate that MARINLINGA is more 
time-efScient compared to Roussopoulos's algorithm, which is well-known for its efficiency. 



m ■ 1 Introduction 
> 

tH- I The line graph I (G) of a graph G is a graph in which every node corresponds to a link of G and two 



nodes are adjacent if and only if their corresponding links are adjacent in G (two links are adjacent if 
they are incident to the same node). The graph G is called the original or root graph of / (G). There 

^D I exist examples of line graphs from social network. Given M clubs and A^ students at an university, 
every student joins two clubs. Each student has different choices (we assume that there are enough 
clubs). We define two networks Gi and G2. The M clubs are the nodes of Gi and two nodes are 

K> ', adjacent if two clubs have the same student as their member. The A^ students are the nodes of G2 
and two nodes are adjacent if two students belong to the same club Clearly, G2 is the line graph of 
Gi. Such pairs (Gi, G2) are common in on-line social networks like Facebook, Twitter and etc., where 
users join the special groups where they share the same interest with others. Computing the line 
graph of a graph and constructing the original graph of a line graph also play an important role in 
link partitioning of communities [6] [1] [5] [10] , bond percolation threshold predictions [18], and it also 
enables us to compare the properties of a random line graph |9] and its original graph. 

The following formula [13] can be used to compute the adjacency matrix of the line graph / (G) of 
a graph G, 

MG) = iR'R),^,-'2i (1) 

where R is the incidence matrix of the undirected graph G. If link j is incident to node i, the entry 
rij of R is 1, otherwise 0. In each column there are exactly two 1-entries. 
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Figure 1: The nine forbidden subgraphs for hne graphs [2]. 



Constructing the original graph is far more complex than computing the line graph. Before con- 
structing the original graph from a given graph, it is important to know whether the graph is a line 
graph. Up till now, the following criteria for a graph to be a line graph exist in the literature: 

• A graph is a line graph if and only if it is possible to find a collection of cliques in the graph, 
partitioning all the links, such that each node belongs to at most two of the cliques (some of the 
cliques can be a single node) and two cliques share at most one node [7]. If the graph is not K^, 
there can be only one partition of this type. 

• A graph is a line graph if and only if it does not have the complete bipartite graph i^i^s as an 
induced subgraph, and if two odd triangleq^ have a common link, the subgraph induced by their 
nodes is the complete graph K4 [T5] . 

• A graph is a line graph if and only if none of the nine forbidden subgraphs (see Figure [TJ is an 
induced subgraph of it [2] . 

• A graph is not a line graph [H] if the smallest eigenvalue of the adjacency matrix ([1]) is smaller 
than —2. 



The complete graph on three nodes K3 is a line graph, which has two different original graphs, 
K^ and Ki^^ (Figure [21] (b)). Except for K^, Whitney's theorem [I7][7] states that, all line graphs 
have only one original graph (isomorphic graphs are considered as the same graph). Based on the 
above criteria and Whitney's theorem, several algorithms for constructing the original graph have 
been proposed [8] [12] [11] [3] . Among those algorithms, Roussopoulos's [12] and Lehot's [8] solutions 
are worth mentioning here. 

Roussopoulos's algorithm starts with choosing an arbitrary link in the input graph and calculating 
the number of triangles containing this link. Depending on this value the starting cell is determined. 
The starting cell is a complete graph K^; if m, = 2 it is a link; if m = 3 a triangle that contains the 



^If every node is adjacent to two or zero nodes of a triangle, it is an even triangle. 



starting link. Having a starting cell of the input graph, the algorithm of Roussopoulos continues to 
find a clique, which is deleted. In addition, in each step the vertices of the clique are labeled by a 
group number. One node in a line graph cannot be assigned to more than two groups (otherwise it 
is not a line graph). The nodes of the original graph are those partitions and all nodes are assigned 
to exactly one partition. In the constructed graph there is a link between two nodes, if the nodes are 
assigned to partitions that have a non-empty intersection. The approach of Roussopoulos is based on 
finding the largest connected components and sequentially the number of triangles that contain this 
link. Theoretically finding the largest connected component is, however, an A^P-complete problem 
|16j . Lehot's solution is based on the characterization of line graphs by Van Rooij and Wilf [8]|15j. 

In this paper, we propose a new algorithm, the Matrix i?elabeling 77V verse LINe Graph Algorithm, 
in short MARINLINGA, that constructs the original graph given the line graph. MARINLINGA 
does not explicitly rely on Whitney's theorem, as all previous companion algorithms, but uses link 
relabeling and endnode recognition in a new way. Via extensive simulation analysis, we have compared 
MARINLINGA with Roussopoulos's algorithm. We demonstrate that MARINLINGA consumes less 
CPU running time. The algorithms are tested on the same machinci and we use the same input line 
graphs for both algorithms. 

2 Link adjacency matrix (LAM) and line graph 

Two nodes of a graph are said to be adjacent if there is a link directly connecting them. The adjacency 
matrix A of a graph contains all information of node adjacency: if node i and node j are adjacent, 
the entry aij = 1, otherwise Uij = 0. Similarly, two links are adjacent if they are incident to the same 
node. 

Definition 1 The link adjacency matrix (LAM) C of a graph G with Nq nodes and Lq links is the 
Lg X Lg symmetric matrix with the entry Cij = 1 if link i and link j of G are adjacent, else Cij = 0. 

The line graph / (G) of the graph G has Nug) nodes and Lug) links, and consequently we have 
Lg = NuG)- According to the definitions of the line graph and the LAM, evidently, the LAM C of G 
is equal to the adjacency matrix Ak^g) of ^ (G), 

C = Ai^G) (2) 

Due to Whitney's theorem and ignoring isomorphisms, for any graph except K^ and i^i^s, one 
can construct the graph exclusively from its LAM. Usually, the (node) adjacency matrix is used to 
represent a graph. Here we use the LAM to specify any graph, except for K-^ and i^i,3. Constructing 
the original graph of a line graph is equivalent to converting a graph representation from the LAM to 
the adjacency matrix. By constructing the original graph directly from the line graph, confusion will 
arise concerning the links in the original graph and the nodes in the line graph. By introducing the 
concept of LAM, we can avoid confusion and facilitate the description of our algorithm MARINLINGA. 



^Processor Intel Core 2 Duo CPU T9600 @ 2.80 GHz and 2.96 GB RAM memory on Java Execution Environment 
JAVA-SE 1.6 and Eclipse IDE (version Galileo 3.5). 
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Figure 2: The four possible configurations in which hnk i is adjacent to hnk j. 



3 Properties of the LAM 

For a simple (undirected, unweighted and without self-loops) graph G{Ng-,Lg) with Nq nodes and 
Lq links, the LAM C has more constraints than the corresponding adjacency matrix A^ besides being 
symmetric and containing only and 1 entries. 

A link i has two endnodes, the left endnode i"*" and the right endnode i~ . Link j also has endnodes 
j+ and j^ . There are four configurations where link i is adjacent to link j, as shown in Figure [2j For 
each single pair of links, the LAM only indicates whether they are adjacent. If they are adjacent, we 
still do not know in which of the four possible configurations this pair of links is adjacent. Fortunately, 
by combining the adjacency relation of 3 or more links, we can determine the configuration of those 
links. 

Definition 2 If m links (m > 2) are adjacent to link i and incident to the same endnode of link i, 
these m links are pairwise adjacent. 

Definition 3 The links, which are adjacent to link i, are defined as the neighboring links of link i. 

Definition 4 The links incident to the left endnode i"*" of a link i are defined as the left-neighboring 
links of i, and the links incident to the right endnode i~ are defined as the right-neighboring links of 
i. 

If we can recognize the link adjacency pattern of a link and its neighboring links, we can specify 
LAM entirely. 

FigureO (a) depicts an example of a link and its neighboring links. The link i has 5 left neighboring 
links at its left endnode i"*", denoted as i+i, • • • , z+s, and 4 right neighboring links at its right endnode 
z~, denoted as z_i, • • • , i-^. The link adjacency pattern of these 10 links is shown in Figure [3] (b). In 
the link adjacency pattern, the labels of the left-neighboring links i+i, • • • ,i+5 are larger than link z, 
and smaller than the right-neighboring links i_i, • • • , i_4. 

Given the configuration of link i and its neighboring links, the corresponding link adjacency pattern 
conforms to the following rules: 

1. the left-neighboring links (such as i+i,--- ,i+5 in the example of Figure [3] (a)) are incident 
to the same endnode i"*", and are said (Definition [2]) to be pairwise adjacent. Similarly, the 
right-neighboring links (such as i_i, • • • ,i_4 in the example of Figure [3] (a)) are also pairwise 
adjacent. This explains the two all-1-triangles (surrounded by the dashed lines) in Figure [3] (b). 
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Figure 3: (a) The configuration of a link i and its neighboring hnks. (b) The corresponding hnk 
adjacency pattern, there is at most one 1-entry in each row/column of the submatrix in yellow. If 
all the entries in green and magenta are 1-entries, the entries of the triangle in white must be also 
1-entries. 



the upper one corresponding to i^ and the second triangle corresponding to pairwise adjacent 
links i-i, ■ ■ ■ ,1-4- 

2. Since there is at most one link between two nodes (multi-links are forbidden), each of the 
left-neighboring links can be adjacent to at most one right neighboring link and vice versa. 
Hence in Figure [3] (b), there exists at most one 1-entry in each row/column of the submatrix 
in yellow. 

We summarize this observation: 

Criterion 5 // the given link adjacency pattern has the following features, it is the link adjacency 
pattern of a link i and its neighboring links (the labels of the left-neighboring links are larger than link 
i, and smaller than the right-neighboring links), 

• All entries of the first row are 1-entries; 

• The triangle bounded by the (n,i+ + l)th column (including the (nj+ + l)th column) is an all-1- 
triangle, where n^+ denotes the number of the left- neighboring links of link i and n^+ > 3; 

• There is at most one 1-entry in each row/column of the submatrix, which is from the 2nd to the 
(nj+ + l)th row and from the {n^+ + 2)th to the {n^+ + n^- + 2)th column, where n^- denotes the 
number of the right-neighboring links; 

• The triangle bounded by the (nj+ + 2)th row (including the (nj+ + 2)th row) is an all-1-triangle. 



Theorem 6 Consider three links i, j and k are pairwise adjacent. If each of the other m links is 
adjacent to all the three links i, j and k, then all the m + 3 links are pairwise adjacent. 
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Figure 4: The smallest forbidden link adjacency pattern. 



Proof. The three links i, j and k are pairwise adjacent and the configuration of i, j and k can be K^ 
or Ki^s, as shown in Figure [2T] (b). If the configuration is K^, other links can be adjacent to at most 
two of i, j and k. However, if the other m links are adjacent to i, j and k, the configuration of i, j and 
k must be i^i^a, and i, j and k have a common endnode. Since each of the m links is adjacent to i, j 
and k, the common endnode of i, j and k must be also an endnode of each of the m links. According 
to Definition [21 all these m + 3 links are pairwise adjacent. ■ 

In Figure [3] (b), links i, i+i and i-^2 are pairwise adjacent, as shown by entries in green. Links i+3, 
z+4 and z+5 are adjacent to i, i+i and i+2, as shown by entries in magenta. By Theorem [6l links i, 
z+i, i+2, i+3, i+A and z+s are pairwise adjacent. 

3.1 The basic forbidden link adjacency patterns in a LAM 

Figure S] (a) depicts the smallest forbidden link adjacency pattern in a LAM. The configuration of 
links i, j and k is a. path on four nodes. Since link i has neighboring links at both of its two endnodes, 
and if link r is adjacent with link i, then link r must be also adjacent with link j or k. Hence, the 
pattern in Figure |4] (a) will not appear in a LAM. 

There are 6 forbidden link adjacency patterns of links i, j, k, r and t, as shown in Figure [H Since 
the number of the left-neighboring links of link i is smaller than 3, we cannot use Criterion [5] to 
prove that the 6 link adjacency patterns are forbidden. However, Figure [5l which exhibits the possible 
configurations of the link adjacency patterns of links i, j, k and r, will facilitate the proof that the 6 
link adjacency patterns in Figure [6] are forbidden. 

The link adjacency pattern of links i, j, k and r in Figure [6] (a), (b) and (c) are the same as the 
link adjacency pattern of links i, j, k and r in Figure[5] (a). There are only two possible configurations 
of this link adjacency pattern. As we can observe in Figure [5] (a), it is impossible to have a new link 
t which is only adjacent with link i, or only adjacent with links i and j, or adjacent with all of i, j, k 
and r. Hence, the patterns in Figure [6] (a), (b) and (c) are forbidden. In the same way, we observe 
that the patterns in Figure [6] (d), (e) and (f) are also forbidden. 

When the number of the left-neighboring links of link i is not smaller than 3 (which implies that 
the number of 1-entries in the first all-1-triangle is not smaller than 6), we can use Criterion [5] to 
determine whether a link adjacency pattern is forbidden. 
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Figure 5: The possible configurations for two link adjacency patterns of 4 links. This figure helps to 
prove that the patterns of 5 links in Figure [6] are forbidden. 
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Figure 6: The forbidden link adjacency patterns of 5 links. 




1 5 10 15 20 25 30 35 40 45 50 
(a) 




1 5 10 15 20 25 30 35 40 45 50 
(b) 



Figure 7: Matrix relabelling on the LAM C of a graph with 50 links. The red dots represent 1-entries. 
(a) Before relabelling; (b) After relabelling. 



4 The matrix relabeling inverse line graph algorithm (MARIN- 
LINGA) 

MARINLINGA is the algorithm that we designed to compute the original graph of a line graph, given 
the adjacency matrix of that line graplu. 

As explained in Section O the adjacency matrix A(g) of ' (G) is equal to the LAM Cq of G. 
Constructing the original graph of a line graph, is equivalent to constructing a graph given the LAM 
of that graph. MARINLINGA only deals with the upper triangle of the given LAM C. 



4.1 Matrix relabeling 

The matrix relabeling algorithm rearranges the LAM C in such a way that the left and right neigh- 
boring links of the first link can be recognized via Theorem [6] and the construction algorithm can work 
efficiently. In each column there are some 1-entries (red dots). If after relabeling the top 1-entries of 
all the columns are connected by a curve, the curve should be nonincreasing. For example, by the 
LAM C of a graph with 50 links in Figure [7] (a), we can only determine which links are adjacent to 
the first link, without any information about which endnode of the first link that the neighboring links 
are incident to. Fortunately, according to Theorem [6l the relabeled LAM C in Figure [7] (b) tells that 
links 2-5 are the left-neighboring links of the first link and links 6-10 are the right-neighboring links. 
Let us first introduce the meaning of swapping the labels of two links in a LAM ClqxLq- The 
entry Cij indicates whether links i and j are adjacent. Swapping the labels of links j and k (j < k) 



■^Although MARINLINGA is designed for connected line graphs, it is also convenient to compute the original graph 
of a disconnected line graph component by component. In the description of MARINLINGA, the connectedness of the 
concerned graph is always assumed. 
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Figure 8: The illustration of swapping the labels of link 5 and 8. The entries in green, magenta and 
yellow ought to be swapped respectively. 



implies that links which are previously adjacent to link j are now adjacent to link k, and links which 
are previously adjacent to link k, are now adjacent to link j, but the adjacency relation between links 
j and k is the same as before, namely the entry Cjk of ClqxLq is unchanged. Hence, swapping the 
labels of links j and k (j < k) means to swap the entries Cij and Cik for i = 1, 2, • • • ,j — 1 (shown in 
the example of Figure [8] in green), the entries Cji and Cj^ for i = j + 1, ■ ■ ■ , k — 1 (in magenta), the 



entries Cji and Cki, i = k + 1,- ■ ■ , Lq — 1, Lq (in yellow). 



Algorithm 1 C <^ SwapLabel{C, j, k) 



1: for i = 1 to j — 1 do 

2: SWap{Cij , Cik) 

3: end for 

4: for i=j + ltok — 1 do 

5: swap{cji,Cik) 

6: end for 

7: for i = A; + 1 to Lq do 

8: swap{cji,Cki) 

9: end for 



Lines 1-3 of the metacode of Algorithm [T] swap the entries Cji and Cik, i = j + 1,- ■ ■ ,k — 1, and 
lines 4-6 swap the entries Cji and Cj^, i = j + 1,- ■ ■ ,k — 1, and lines 7-9 swap the entries Cji and c^i, 
i = k + 1,- ■ ■ , Lg — 1, Lq. The code swap {cij,Cik) of line 2 is equivalent to the codes: t = Cij; Cij = Cik] 

Next, we will explain the matrix relabeling algorithm. We will first give an example showing how 
the matrix relabeling algorithm relabels the LAM C in Figure [7] (a) into the matrix in Figure [7] (b). 
In the first row of the matrix in Figure [7] (a) there are 9 1-entries in total. There are 6 0-entries 
from ci,2 to ci,io and 6 1-entries from ci,ii to ci,5o: ci,3 = ci,5 = ci,6 = ci,8 = ci^g = ci,io = and 
ci^i3 = ci^is = ci^is = ci^ig = ci^24 = ci^4o = 1. We swap the labels of links 3 and 13, links 5 and 




1 5 10 15 20 25 30 35 40 45 50 



Figure 9: The LAM C after the relabehng of the first row. 



15, hnks 6 and 18, links 8 and 19, hnks 9 and 24, hnks 10 and 40 by Algorithm [2] and the LAM C is 
shown in Figure [H In the second row, there are 3 1-entries from C2,3 to C2,io- There are 2 0-entries 
from C2,3 to C2,5 and 2 1-entries from C2,6 to C2,io: C2,4 = C2,5 = and C2,6 = C2,9 = 1- We swap the 
labels of links 4 and 6, links 5 and 9. By similar operations, we relabel the LAM C into the order 
shown in Figure [7] (b) . 

Now we give the general description of the matrix relabeling algorithm. In the fcth row of C, Lines 
1-7 of Algorithm [2] store the value of i in X when the entry c^i is 0, i = ti + 1, • • • ,a + u. Lines 8-14 

h 

store the value of i in y when the entry c^j is 1, i = a + u + 1, • • • , 6. If a = ^ c^j, X and Y have 

i=u+l 

the same number of elements. Lines 15-17 swap the labels of Xi and 1^, where Xi and Yi are the iih 

element of X and Y respectively. For instance in Figure [8] (b), if we take n = 2, A: = 2, 6 = 10 and 

10 r 1^ r 

a = X;c2i = 5, by Algorithm^ X = 5 7 ,y= 8 10 



the labels of links 5 and 8, 7 and 



10 are swapped respectively. 

Lines 1-2 of Algorithm [3] make the neighboring links of link 1 have the smaller labels than the 
other links. By lines 3-4, the labels of the links which are adjacent to both link 1 and 2 are smaller 
than those of the remaining links. Further, lines 5-6 let the labels of the links which are adjacent to 
all of links 1, 2 and 3 are smaller than those of the remaining links. Lines 7-14 make that the labels 
of the links which are adjacent to link i but not adjacent to links 1, ■ ■ ■ ,i — 1, are smaller than the 
labels of the links which are not adjacent to link 1, • • • , i, for i = 2, • • • , Lq- Figure [7] and [10] show 
examples of C before and after matrix relabeling. 

Lg si+l S2+2 

Let si = ^cii, S2 = X] C2j and S3 = ^ csj. After relabeling by Algorithm [3l the given LAM C 

1=2 i=3 i=4 

satisfies: 



10 



Algorithm 2 C <^ GroupLabelSwapping{C, u, k, a, b) 



1: m <^ 

2: for i = u+l to a + n do 

3: if Cki = then 

4: m -^ ra+1 

o. y\^Yfi "^ — '' 

6: end if 

7: end for 

8: m ^ 

9: for i = a + n+lto6do 

10: if Cki = 1 then 

11: m -^ ra+1 

13: end if 

14: end for 

15: for i = 1 to ra do 

16: C <^ SwapLabel{C, Xi,Yi] 

17: end for 



Algorithm 3 (C, si,S2,S3) ^ MatrixRelabeling{C) 



= the sum of cij, where i = 2 to Lq 
= GroupLabelSwapping{C, 1, 1, si, Lq) 
= the sum of C2i, where i = 3 to si + 1 
= GroupLabelSwapping{C, 2, 2, S2, si + 1) 
= the sum of c^i, where i = 4 to S2 + 2 
= GroupLabelSwappingi^C, 3, 3, S3, S2 + 2) 
si + 1 

:2 

9: while s < Lq and A; < Lq do 

10: s <^ the sum of c^i, where i = s + 1 to Lq 
11: C <^ GroupLabelSwapping{C, s, k, s, Lq) 
12: k ^ k + 1 

13: S^S + S 

14: end while 
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(a) 




1 100 200 300 400 508 
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1 100 200 300 400 462 




1 100 200 300 400 496 



Figure 10: The relabeled C of four ER random graphs G {N,p): (a) N = 350, p = ^^f^; (b) A^ = 200, 
p = ^-^f^; (c) N = 100, p = ^^^^; (d) N = 32,p = l, where p = ^^^ is the threshold probability 
for the connectivity of the graph. 
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Figure 11: The LAM (a) relabeled by Algorithm [3] and its corresponding graph (b). 

• For i = 2,- ■ ■ , si + 1, cij = 1; and for i = si + 2, • • • , Lq, cu = 0. 

• For z = 3, • • • , S2 + 2, C2i = 1 if S2 > 1; and for i = S2 + 3, • ' ' ) •5i + 1, C2j = if si > S2 + 2. 

• For z = 4, • • • , S3 + 3, C3i = 1 if S3 > 1; and for z = S3 + 4, • • • , S2 + 2, C3j = if S2 > S3 + 2. 

• If link j {j > si+i) is adjacent to link i but not adjacent to links 1, 2, • • • ,i — 1 (i > 2), and link 
k (k > si+i) is not adjacent to all of links 1, 2, • • • , i {i >2), then j < k. 

If 53 > 1 (which implies that S2 > 2 and si > 3), according to Theorem [H links 2, 3, • • • , S3 + 3 are 
the left-neighboring links of links 1 and the links S3 + 4, • • • , si + 1 are the right-neighboring links of 
link 1, as illustrated in the example of Figure [TT] where si = 9 and S3 = 3. 

4.2 Construction algorithm 

The construction algorithm converts the relabeled C into the matrix £'2xLq, where the entries en and 
e2i denotes the two endnodes of link i. During the process of the construction, the zero entries of 
E2xLg mean that the endnodes have not been determined yet. 

Section 14.2.11 will first show an example of graph construction, and section 14.2.21 and 14.2.31 will 
describe the general construction algorithm. 



4.2.1 An example of graph construction from C 

From the given LAM C in Figure [7] (b), we deduce that the graph has 50 links. Based on the LAM 
C, we will determine the endnodes of the 50 links. The construction consists of the following steps: 

1. Let nodes 1 and 2 be the endnodes of link 1. According to Theorem [6l node 1 is also the endnode 
of links 2-5 and node 2 is also the endnode of links 6-10, as shown in Figure [12] (a) and equation 
([3]) below, where the numbers above the matrix are the link numbers. 
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Let node 3 be the other endnode of hnk 2. The 2nd row of the LAM C shows that hnks 11-14 
are adjacent to hnk 2. Hence, node 3 is also the endnode of hnks 11-14, as shown in Figure [12] 
(b) and equation @. 
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Similarly, let node 4 be the endnode of link 3, 6 and 15-18 as shown in Figure [12] (c) and equation 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 




1 


1 


1 


1 


1 


2 


2 


2 


2 


2 


3 


3 


3 


3 


2 


3 


4 








4 



























E 



15 


16 


17 


18 


19 




50 


4 


4 


4 


4 
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and let node 5 be the endnode of link 4, 8 and 19 as shown in Figure [T^ (d) and equation (|H|), 
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(6) 



and let node 6 be the endnode of link 5, 16 and 20-23 as shown in Figure [12] (e) and equation 

dZD- 
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(7) 



Then compute the LAM of the constructed part of the graph as shown in Figure [13] The red 
dots are 1-entries which are from the given LAM in Figure [7] (b) . The green dots are 1-entries 
which are determined by the red 1-entries. If the corresponding entries in the given matrix are 
not 1, then the matrix is not a LAM. 

2. In the second step, we scan rows 6 to 10 of the LAM, since links 6 to 10 are incident to the same 
endnode. Let node 7 be the endnode of link 7, 21 and 24-25, and let node 8 be the endnode of 
link 9 and 20, and let node 9 be the endnode of link 10, 14 and 26-27, as shown in Equation ([8|) 
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Figure 12: The example of construction. The initiahzation is done in (a). Both or one of the two 
endnodes of Hnks 1-23 are determined. 
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Figure 13: The LAM of the constructed part (Hnks 1-23) of graph are computed. The green 1-entries 
are determined by the red 1-entries. 
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Figure 14: The example of construction. Both or one of the two endnodes of hnks 1-27 are determined. 



and Figure [Ml 
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3. Similarly, let node 10 be the endnode of link 11, 19 and 28-30, and let node 11 be the endnode 
of link 12, 18 and 31-35, and let node 12 be the endnode of link 13 and 36, as shown in Equation 
dSI and Figure [161(a). 
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Figure 15: The LAM of the constructed part (links 1-27) of graph are computed. The green 1-entries 
are determined by the red 1-entries. 




Figure 16: The example of construction. Both or one of the two endnodes of links 1-36 are determined. 



17 




10 15 20 25 30 35 



Figure 17: The LAM of the constructed part (hnks 1-36) of graph are computed. The green 1-entries 
are determined by the red 1-entries. 



4. Constructing in this way, the two endnodes of ah the hnks are eventuaUy determined, as shown 
in Equation (|1U|) and Figure [TH] (a) . The final structure of the matrix E exhibits the hnk hst 
of the original graph G which consists of 30 nodes and 50 links. For example, link 36 connects 
node 12 and node 15 in G. The matrix E is readily transformed into the adjacency matrix of G. 
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4.2.2 Initialization (The recognition of the endnodes of the first link and its neighboring 
links) 

When S3 > 1, Theorem [6] implies that S2 > 2, si > 3 and links 2, 3, • • • , S3 + 3 are incident to the 
left endnode of link 1 and links S3 + 4, • • • , si + 1 are incident to the right endnode of link 1. Therefore, 
line 1-2 of Algorithm [H initialize E by £. The numbers above the matrix £ in (llip are the column 
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Figure 18: The example of construction. The two endnodes of ah hnks are determined. 
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Figure 19: The LAM of the constructed graph is computed. The green 1-entries are determined by 
the red 1-entries. 
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Algorithm 4 E2xLg -^ Initialization{C, 81,82,83) 



if S3 > 1 then 

E^£ 
else if si = 1 then 

E ^£1 
else if si = 2 then 

6: ^2xLg -^ Initialization2{C,82,83) 

7: else if si = 3 then 

8: E2xLg <= I'n'itialization3{C,82,83) 

9: else if si > 4 then 

10: E2xLg -^ Initialization4:{C,8i,S2,S3) 
11: end if 



numbers, which indicate the Unk numbers, and £ has the following structure, 

1 2 ■■■ S3+3 S3+4 ••■ si+1 si+2 ••• Lg 
I 1 I 1 I ••• 
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(11) 



When S3 = 0, Theorem [6] cannot be used. However, the limited number of cases of S3 = enables 
us to accomplish the initialization with the detailed analysis in the Appendix Rl 

4.2.3 The recognition of the endnodes of the whole graph 

Lines 1-2 of Algorithm [5] relabel the given LAM C and determine the initial state. In the initial state, 
link 1 is always incident to node 1 and 2. Some of the neighboring links of link 1 are incident to node 
1, and the other neighboring links are incident to node 2. The second endnodes of the neighboring 
links of link 1 have not decided yet in the initial state. 

Line 3 initiates the number of nodes Nq to 2. The two endnodes of link 1 are already determined. 
Starting with link 2 until link Lq (line 4) , the number of nodes Nq increases by 1 (line 6) if the second 
endnode of link i is not determined (line 5). Let the second endnode of link i be Nq (line 7). When 
link i is adjacent to link j, j = i + 1,- ■ ■ , Lq (lines 8-9), let the first endnode of link j be Nq (line 11) 
if the first endnode of link j is not determined (line 10) . If the first endnode of link j is determined 
but the second endnode is not determined and links i and j do not share the first endnode (line 12), 
let the second endnode of link j be Nq (line 13). 

4.3 Worst case complexity of MARINLINGA 

Algorithm [1] has a complexity of O (Lg)- The complexity of Algorithm [3] can be computed as follows. 
Line 1 has a complexity of 0{Lg)- In the worst case, the function of line 2, Algorithm [2] has a 
complexity of O {Lq^ , if m in line 15 of Algorithm [2] is proportional to Lq- The worst case complexity 
of lines 3-6 is also O {Lq). Hence, lines 1-6 have a complexity of O [L'qJ. Neglect O (1) operations of 
lines 7-8. The times that lines 9-14 are executed is stored in k. If k is proportional to Lq, rri in line 
15 of Algorithm [2] niust be bounded by a constant, then the complexity of line 11 is O (Lg)- If k is 
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Algorithm 5 E^axLc <= MARINLINGA{C) 



1: 


(C, Si, S2, S3) <^ Matrix Relabeling{C) 


2: 


E2xLg ^ Initialization{C, si, S2, S3) 


3: 


iV^2 


4: 


for i = 2 to Lg do 


5: 


if e2j = then 


6: 


N ^N + 1 


7: 


e2i^N 


8: 


for j = i + 1 to Lg do 


9: 


if Cij = 1 then 


10: 


if eij = then 


11: 


eij 4= N 


12: 


else if e2j = and en ^ eij then 


13: 


e2j 4= N 


14: 


end if 


15: 


end if 


16: 


end for 


17: 


end if 


18: 


end for 



bounded, the complexity of line 11 will be O {Lq) . Therefore, lines 9-14 have a worst case complexity 
of O {Lq). Hence, the complexity of Algorithm [3] is O {Lq). 

Algorithm El [7] and [5] have a worst case complexity of O (1), hence the complexity of Algorithm 
m is also O (1). Lines 4-18 of the main Algorithm [5] have a worst case complexity of O {Lq^. In 
summary, the worst case complexity of the MARINLINGA is O {Lq^ . Since the number of links of 
the original graph G and the number of nodes of the line graph I (G) are equal, Lq = Nuq\, the worst 



case complexity of the MARINLINGA is written as O ( N? 



1(G) 

5 Comparison with Roussopoulos's algorithm 

We use the same input line graphs for both MARINLINGA and Roussopoulos's algorithm. We start 
with line graphs constructed from Erdos-Renyi random graphs Gp {Nq) [1]. We calculate the prob- 
ability density function of the difference AT between the running time of Roussopoulos's algorithm 
(TroussopouIos) and MARINLINGA (Tmarinlinga) 

AT = TroussopouIos - ^MARINLINGA (12) 

We randomly create 1000 different line graphs based on the class of Erdos-Renyi random graphs 
Gp (No) for each number of nodes Nq = {10, 20, 30, 50} and link density p = {0.1, 0.2, ..., 0.9}. The 
probability density functions of the time difference /at {x) for each class of line graphs of Gp (No) 
are shown in Figure [20j The values of the probability density function are nearly always positive 
which means in practice that MARINLINGA needs less time for the execution than Roussopoulos's 
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Figure 20: PDFs of the AT for N = 10, 20, 30 and 50 and p = {0.1, 0.2, . . . , 0.9} 
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Ng 
P 


10 


20 


30 


50 


.1 


0.0462 


0.7839 


3.6660 


42.8153 


.2 


0.1682 


4.3012 


26.1086 


385.3208 


.3 


0.6281 


11.9492 


49.826 


1582.5827 


.4 


1.3921 


26.4580 


209.3056 


4802.8110 


.5 


2.4269 


51.8393 


422.8641 


7060.6021 


.6 


3.5808 


91.8258 


719.6185 


10978.1532 


.7 


4.9181 


145.5864 


1165.1978 


15942.7395 


.8 


6.8952 


217.2041 


1768.9901 


23393.5148 


.9 


9.2986 


317.4966 


2501.4799 


24723.4739 



Table 1: Expectations of the time difference (//sec) 

algorithm. 

We calculate the expectation according to [l3] and the experimental results for all of the mentioned 
cases AT,- for i = 1, 2, . . . , 1000. 



E[AT] = ^k Pi[AT 



k] 



V^lOOO A.rTi 

1000 



(13) 



The results in milliseconds are given in the Tabled) 

Additionally, we calculate the probability that MARINLINGA is slower than Roussopoulos's algo- 
rithm: Pr[AT < 0] for each Ng and p. The simulation shows that Pr[AT < 0] > only for Ng = 10 
and p < 0.2, in which the graphs are mostly disconnected. When the graph is disconnected, MARIN- 
LINGA needs extra time to partition the graphs into connected components (see footnote 3). For 
Ng = 10 and p = 0.1, Pr[Ar < 0] = 0.33 and for Ng = 10 and p = 0.2, Pr[AT < 0] = 0.01. For all 
the other cases 



Pr[Ar < 01 < 0.001 



(14) 



which means that MARINLINGA is generally more efficient than Roussopoulos's algorithm. The 
algorithm to find the maximal connected common subgraphs in graphs is frequently used in the 
Roussopoulos's algorithm. This algorithm requires a high running time, because the problem of 
finding the maximal connected common subgraphs is A^P-complete |16] . The dependence on this 
A^P-complete algorithm is most significant weakness of Roussopoulos's algorithm. 



6 Conclusion 

We have presented a new algorithm MARINLINGA for reverse line graph construction. By introducing 
the concept of LAM, we transformed the problem of reverse line graph construction into the problem 
of constructing a graph from the LAM. MARINLINGA consists of two sub-algorithms: the matrix 
relabeling algorithm and the construction algorithm. The matrix relabeling algorithm preprocesses 
the LAM into the special order by which we can determine the neighboring links of the first link and 
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the endnodes of the first hnk incident to the neighboring links. The construction algorithm makes the 
first two nodes be the endnodes of the first link by default, and thereafter, determines the endnodes of 
the remaining links. MARINLINGA has a worst case complexity of 0{N^,q-.), where Niiq\ denotes the 
number of nodes of the line graph. We have demonstrated that MARINLINGA is more time-efficient 
compared to Roussopoulos's algorithm for connected line graphs. The complexity of Roussopoulos's 
algorithm mentioned in [12] is 0{Ni(^q^ +-^i(G))) where Nk^q^ and -^^(g) are number of nodes and links of 
the line graph. Since ^^(g) = O \Nf,Q^ J in worst case, the complexity of Roussopoulos's algorithm is 
also 0{Nf,^C) in worst case. However, this analysis neglects the computational time of a sub-algorithm 
that determines the maximal connected common subgraph in each iteration. Finding a maximally 
connected common subgraph is an iVP-complete problem [16J, implying that Roussopoulos's algorithm 
is, in fact, not polynomial in worst case. 
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A The initialization of the construction algorithm when S3 = 

Theorem [6] cannot be used when S3 = 0. Since there exists limited number of cases of 53 = 0, we can 
still accomplish the initialization. 

A.l When si = 1 

Link 1 has only one right neighboring link: link 2. Link 1 does not have left neighboring links. The 
initial state of E is £i. Lines 3-4 of Algorithm [J] initialize E hy £i. 



£i 



1 2 ••• 

2 ••• 



(15) 



A. 2 When si = 2 



There are different adjacency patterns. The submatrix of C in Figure [21] (a) implies that, links 2 and 
3 are adjacent to link 1, and link 2 is not not adjacent to link 3. Links 2 and 3 must be incident 
to two different endnodes of link 1. The pattern in Figure [2T] (b) has two possible configurations Kj, 
and i^i,3. If si = 2 and S2 = 0, the initial state is £2,a^ as shown in lines 1-2 of Algorithm [UJ When 
si = 2 and S2 = 1, because the graph is connected, either C2,4 = 03^4 = 1 or 02,4 = 0,03^4 = 1 or 
C2,4 = 1,C3^4 = 0. If C2,4 = 03^4 = 1, the initial state is £"2,6.1, which is K^, otherwise the initial state is 
£2,b.2, which is Ki^3, as shown in lines 8-12 of Algorithm [6l 
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Figure 21: The adjacency patterns of link 1 and its neighboring Unks when si = 2. The graphs on the 
right are the possible configurations correspondingly. 
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Algorithm 6 E2xLg -^ Initialization2{C , S2) 



1: if S2 = then 

2: E^ S2,a 

3: else 

4: if C2,4 = 1 and 03^4 = 1 then 

5: E -^ £"2,6.2 

6: else 

7: E^ ^2,6.1 

8: end if 

9: end if 



A.3 When si = 3 

There are two recognizable adjacency patterns as described in Figure [22] (b) , and (c). Taking pattern 
(c) as an example, links 1, 2 and 3 are pairwise adjacent, then the configuration of them is K3 or 
i^i,3, as shown in Figure [21] (b). Link 4 is also adjacent to link 1, but not adjacent to links 2 and 3, 
suggesting that the configuration of links 1, 2 and 3 must be K3, and link 4 is incident to the other 
endnode of link 1. Figure [22] (a) depicts the smallest forbidden link adjacency pattern in a LAM. The 
adjacency relation of links 1, 2 and 3 is recognizable, and the configuration is a path on four nodes, 
as shown in Figure [21] (a) . Link 4 is adjacent to link 1, then link 4 must be also adjacent to links 2 or 
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Figure 22: The adjacency patterns of link 1 and its neighboring hnks when si = 3. Pattern (a) is 
forbidden, and patterns (b), (c) and (f) correspond to only one configuration respectively. Patterns 
(d) and (e) both have two possible configurations. 



3. Hence the pattern is forbidden. If si = 3 and S2 = 0, the initial state is £'3_b (lines 1-2 of Algorithm 
[7]). If si = 3, S2 = 1 and 03^4 = 0, the initial state is £^^c (lines 3-5 of Algorithm [7]) . When si = 3, 
S2 = 1 and C3^4 = 1, due to the connectivity of the concerned graph, either 02,5 = 03^5 = 04^5 = 1 or 

C2,5 = C3,5 = 1, C4,5 = Or C2,5 = C3,5 = 0, €4,5 = 1 Or C2,5 = 1, C3,5 = €4,5 = Or C2,5 = 0, 03,5 = €4,5 = 1. 

If C2,5 7^ C3^5, the initial state is £"3,^.2, else if 02,5 = 03^5 7^ 04^5, the initial state is 1^3,^.15 else if 
C2,5 = C3,5 = C4^5 = 1, we need to look further at the relation of C2,6 and c^^^: if C2,6 7^ c^fi, the initial 
state is £"3,^.2, else the initial state is £"3,^.1 (lines 11-15 of Algorithm [T]) . If there are only 5 links in 
total and 02,5 = 03,5 = 04,5 = 1, one can choose any of £"3,^.1 and £^4.2 as the initial state, and get 
isomorphic configurations. If si = 3, S2 = 2 and S3 = 0, the same method is employed (lines 21-26 of 
Algorithm [7|) . 
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A. 4 When si > 4 
A.4.1 When S2 > 3 

The configuration is unique. The initial state of E is £4^. 
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Algorithm 7 E2xLg -^ Initialization3{C,S2,s^) 



1: if S2 = then 

2: E^ S3,b 

3: else if S2 = 1 and 03^4 = then 

4: E<^ £3,0 

5: else if S2 = 1 and 03^4 = 1 then 

6: if C2,5 / C3,5 or (c2,5 = £3,5 and 02,5 = 04,5 and Lq = 5) then 

7: ii^ <^ ^3,d.2 

8: else if 02,5 = 03^5 and 02,5 7^ 04^5 then 

9: E <^ e^^d.i 

10: else if 02,5 = 03^5 and 02,5 = 04^5 and C2,6 = £3,6 then 

11: E <= e^^d.i 

12: else if C2,5 = £3,5 and 02,5 = £4,5 and C2,6 7^ £3,6 then 

13: £; <= &i,d.2 

14: end if 

15: else if S2 = 2 and S3 = then 

16: if C2,5 / C3,5 or (c2,5 = C3,5 and C2,5 = C4,5 and Lg = 5) then 

17: E 4= ^3,e.2 

18: else if 02,5 = 03^5 and 02,5 7^ 04^5 then 

19: E 4= Ss^e.l 

20: else if 02,5 = 03^5 and 02,5 = 04^5 and C2,6 = £3,6 then 

21: E <= Ss^e.l 

22: else if C2,5 = C3^5 and 02,5 = 04^5 and C2,6 7^ £3^6 then 

23: E 4= £"3,6.2 

24: end if 



25: end if 
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Figure 23: The adjacency pattern and the corresponding configuration when S3 = 0, S2 > 3 and 

si >4. 
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A.4.2 When S2 < 2 

There are 13 forbidden patterns, as shown m Figure [2^ where the hnks with labels larger than 5 are 
not displayed. The pattern in Figure [22] (a) is forbidden, hence the 4 patterns in Figure [Ml (a.l) are 
also forbidden, where x can be 1 or 0. The pattern of links 1 — 4 in Figure [Ml (a. 2-3) is the same as 
the pattern in Figure [22] (b), which has a specific configuration. In Figure [Ml (a-2), link 5 is adjacent 
to link 1 but not 2, then link 5 must be adjacent to link 3, which is not true, hence the 2 patterns 
in Figure [M] (a. 2) are forbidden. In Figure [Ml (a. 3), link 5 is adjacent to link 1 and 3, then link 5 
must be adjacent to link 4, which is not true, hence the pattern in Figure [Ml (a. 3) is also forbidden. 
Similarly, based on the patterns in Figure [22} we can conclude that patterns in Figure [Ml (b.l), (b.3), 
(c.l), (c.3), (d.l) and (d.4) are also forbidden. Based on the values of entries S2, 03^4, 03^5 and 04^5, 
Algorithm [8| decides the initial state of E. 
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Algorithm 8 E2xLg -^ I'nitialization4{C, si, S2, s^) 



1: if S2 > 3 then 

2: E <^S4 

3: else if S2 = or (s2 = 1 and 03^4 = 1 and 03^5 = 1 and 04^5 = 1) then 

4: E^ £■4,^.4 

5: else if S2 = 1 and 03^4 = and 04^5 = 1 then 
6: E <^ £'i,b.2 

7: else if S2 = 1 and 03^4 = 1 and 03^5 = and 04^5 = 1 then 
8: E ^ £'4,c.2 

9: else if S2 = 2 and 04^5 = 1 then 
10: E <= £i^d.2 

11: else if S2 = 2 and 04^5 = then 

12: E <^ £'4,d.3 

13: end if 
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Figure 24: The adjacency patterns of link 1 and its neighboring hnks when si = 4. There are 16 
forbidden patterns. The other 12 possible patterns correspond to only one configuration respectively. 
The entry x can be 1 or 0. 
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